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Geomagnetic field lines are traced from points on a ballistic 
trajectory to the points where they cross a given altitude shell. 
The intersects thus obtained are called ’’trajectory determined 
intersects." Three parametric angles identify and describe their 
positions: the look-angles ^ande , azimuth and elevation respec- 
tively, which relate and intersect to a fixed location on the earth; 
and the line-angle a, which for an observer stationed at this fixed 
location is a measure of the divergence of the line -of- sight from 
the field line. The locus of all trajectory-determined intersects 
is the "intersect trace" at the specified altitude level. 

The intersect trace of a sample rocket trajectory is obtained 
for the 100 kilometer level. The variations of the parametric 
angles with time and altitude are calculated for 19 ground loca- 
tions distributed symmetrically around the launch site. For the 
specific flight path considered, there exists an area on the earth 
where optimum parametric angle conditions prevail, i.e., mini- 
mum changes ini p and e, and smallest values of a. 

Finally, the observer-intersect distance is calculated for all 
stations. 
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PARAMETRIC ANGLES OF TRAJECTORY-DETERMINED 
FIELD LINE INTERSECTS 


by 

E. G. Stassinopoulos 
Goddard Space Flight Center 


INTRODUCTION 


An interesting and promising development in magnetospheric research is the introduction of 
controlled geophysical experiments. 

The first proposal for such an experiment was made by Hess (1965), who suggested that a 
stream of electrons generated by a space -borne accelerator, could produce artificial auroras. 
These auroras would be powerful tools in understanding magnetic field line geometry, conjugate 
point locations, large scale magnetic and electric fields, auroral spectroscopy, plasma instabili- 
ties, and VLF generation. To evaluate the feasibility of the project, Hess suggested in a prelim- 
inary outline of a testing plan (private communication, 1967), that the initial apparatus be flown on 
an Aerobee 350 from Wallops Island. 

In this report we investigate topological aspects of the experiment and appraise the position 
and orientation problems associated with visual or instrumental tracking of artificial auroral 
spots. We have chosen for study a nominal Aerobee 350 trajectory with a launch azimuth of 150 
degrees, launch elevation of 87 degrees, and effective payload of 611 pounds. 

To observe artifically generated auroras, especially with instruments of relatively small view 
angles, it is desirable to have advance knowledge of the approximate location where the auroral 
spots are expected to appear. If the particle source happens to move along a ballistic trajectory, 
as in the proposed experiment, the observational difficulties are compounded because to nearby 
observers the auroral path may sometimes look erratic. In fact, it will fluctuate strongly in 
azimuth and/or elevation if a station is within a "critical" range around the launch site. Rapid 
displacements of the spots in the field of view of the instruments are then to be anticipated. 


Sudden and swift variations of the look angles for close -by observers could seriously impair 
the success of a costly experiment. A search-and-aim operation in real time could prove futile 
without some detailed prediction, owing to the short duration of actual beam-on intervals. This 
obstacle can be overcome if the operator anticipates these changes and knows beforehand in what 
direction they will move. He could, in fact, practice the motions in a sort of dry run, long before 
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the real launch. For this purpose a program has been designed which can approximately simulate 
such an experiment. Sample calculations were made with a theoretical flight path. Lines of force 
passing through selected points on this path were traced in the direction of the field. Their inter- 
sects with a set of successive altitude levels or "shells", ten kilometers apart, called "trajectory- 
determined intersects”, were determined. 

Three parametric angles identify and describe the position of those intersects: the look- 
angles i/> and 6 , azimuth and elevation respectively, which relate an intersect to a fixed location on 
the earth; and the line -angle a, which for an observer stationed at this fixed location is a measure 
of the divergence of the line-of-sight from the field line. 

The "intersect trace" or locus of all such intersects for the 100-km shell was obtained. 
Parametric angles for 19 strategically placed stations were computed. The results were plotted 
versus time and the "critical" range was established. A region of optimum-parametric-angle- 
conditions emerged. 

PARAMETERS OF PROBLEM 

Trajectory, L-band, Time Limitations 

The trajectory used for the computations is within the launch capability of the Wallops Island 
facilities; it pertains to an Aerobee 350 rocket with an effective payload of 611 pounds launched at 
150° azimuth and 87"’ elevation. Owing to the low latitude of the launch site and the southeasterly 
direction of the impact zone, it crosses relatively few ”L-shells". The L-band traversed by the 
flight path is shown in Figure 1; it has a mean value of about 2.56 earth radii. 

Field lines were traced only for points on the trajectory above 150 km. This limited the study 
to the time interval bounded by t' = 110 and t' = 460 seconds, where the prime indicates "after 
liftoff". The evaluation of the parametric angles was further restricted to 200 < t' < 400 seconds 
because during the initial phase of the flight, that is for t' < 200 seconds, unsettled conditions 
prevail. 

Specifically, until its burn-out at about t' = 55 seconds, the rocket is affected by perturbations 
which may cause the actual trajectory to deviate substantially from that originally planned. An 
accurate prediction of the true path to be followed by the rocket cannot be made on the basis of 
initial conditions, such as launch azimuth or elevation, payload, and thrust. Maneuvers and 
activities, required to stabilize the vehicle, orient the apparatus to the field direction, and deploy 
the ion collecting screen usually would be scheduled during the 145 seconds following burn-out. 
Furthermore, any other task or function that may be necessary prior to the activation of the ex- 
periment would be performed then. 

During the terminal phase of the flight, when t' > 400 seconds, the increasing speed of the 
descending vehicle hampers effective tracking of the spots with instruments of narrow view angles. 
Therefore, with regard to parametric angles, the best time for conducting experiments lies between 
200 and 400 seconds after liftoff. 
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Stations 


Figure 2 depicts the location of the nineteen stations selected for investigation; it also shows 
the launch site, the horizontal projection of trajectory points 10 seconds apart, and the correspond- 
ing intersects at 100 km altitude; furthermore, the position of the peak altitude point is indicated 
and the beginning of the "prime" time interval for experiments is marked. 

At first, eleven reference stations were placed in a circle around the peak altitude intersect. 
This arrangement, yielding more information on the look-angle-to-position relationship than any 
other pattern, also saves valuable computer time. Later, eight more stations were added in a 
linear formation to investigate the parametric angles along a line parallel to the main segment of 
the intersect trace. 

METHOD 

Field Line Tracing and Altitude Intersects 

A method developed by the author for geomagnetic field line tracing was presented in a pre- 
vious paper (Roederer et al., 1966). A brief review of this method is given here in Appendix A. 

The procedure used in this report is essentially the same, with only minor modifications to adapt 
it to the present work.* 

The magnetic field strength and field vector components are calculated with Mcllwain's new 
MAGNET subroutine (part of the INVAR subroutine), employing the 99-term geomagnetic field 
model by Hendricks and Cain (1966) for the epoch 1960.0, updated to 1965.0. f The L-parameter 
is computed with Hassit and Mcllwain's new (1967) INVAR subroutine. 

The field line integration stepsize, a constant of the tracing process, was set to 2 kilometers; 
tests with larger elements of arc had produced disproportionate errors in the position of the inter- 
sects, up to ±1 degree in latitude and/or longitude, while still expending almost the same computer 
time. 

Lines were traced in the direction of the field only. 

The Look-Azimuth ^ 

The azimuth of a point in space C relative to a position a on a sphere is the angle between the 
meridian plane through a and the plane S(a, b, 0). measured clockwise from the northern direction 


•The following changes were introduced into the line tracing procedure: The direction of the tracing was not reversed, since no con- 
sulate r.rrArA- thr «::arstit:es ' 'arclencth” ? “din” and “declination” were of no concern and were not computed; no 

search for the existence and location of a minimum-B point was made; the L computing part was retained without the value test; the 


section that calculates the intersects was extended to permit their calculation at several consecutive altitude levels during each 
tracing operation. 

tThis model is best known as the GSFC 9/1965. It is available from the National ^pace Sciences Data Center at Goddard Space Flight 
Center, Greenbeit, Maryland 20771. 
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of the reference meridian (Figure 3a). The angle p in Figure 3b is the azimuth of c with respect 
to a. From spherical trigonometry: 


( cos 8 c - cos 0 cos p\ 

s c in P sin p J = Up, 0. e c ) • (l) 

For the geoid, when a and b are related to the known vectors r and R c (Figure 1), 8 and 6 c are 
constants and Up, 0 C ) is a function of p only. The cosine of p can be readily found from: 


cos p 


R • R 

c 

RR 


i 8 sin # c cos “ <^ c j + cos 0 cos 0 c 


( 2 ) 


Thus, p is a function of the four position angles of a and b. The azimuth of a field line intersect is 
easily determined from Equation 2. 

The Look-Elevation € 

The elevation e of a point in the sky with respect to an observer on the surface of the earth is 
the angle formed by the horizontal plane K and the line from the observer 0 to the point c as in 
Figure 4a. In a geocentric spherical coordinate system it is convenient to evaluate the elevation 
vectorially, in terms of the variables R and R c and the angle p between them. 

Let Q be some vector along the projection of c in the horizontal plane K, pointing away from the 
observer (Figure 4a). Obviously, in unit vector notation, 

sin e “ C • R . (3) 

cos e = C • Q . (4) 

We see from Figure 4b that 

C = 

and from Figure 4c and 4d that 


|ch R - R - 


RR 


(R x R c ) x R 


(5) 


(61 
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These relations give: 


£ - sin ‘ |t=t R c cos p - R 


= COS " 1 R c sin p 


(?) 


The Line-Angle a 

In the southern magnetic hemisphere, the angle formed by a descending field line B and the 
line of sight c of an observer will be called the "line angle" and will be represented by a (Figure 
5a). Obviously 


„ C • B 

COS p - = - COS CL . 


(8) 


Since 


C = R c -R = f(R Cx -R x ) + j( R c y - R y) + k( R c 2 -R 2 ) 


(9) 


the cartesian components are 


C = R - R = R sin 6* cosdj - R sin 8 cos <j> , 

X C X c c 


C - R - R R sinfr sin<£ - R sin 6 sin 4> , 

v c y c c T c 

y y 


C z = R c - R z - R c cos6' c - R cos 0 . 




( 10 ) 


Since a unit vector conversion in Equation 9 leads to 


C r = C x sin 6 c cos 4> c + C y sin & c sin <t> c + C z cos 6 ^ . 


Cg = C x cos ^ c cos 4> c + C y cos0 c sin cf- c - C z sin£? c 


C x sin i c + C y cos 4> c 


(ID 


J 
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the spherical components are 


c 

r 

= c 

■ ? 

R c - R cos 4 > , 




c e 

- c 

■ 6 = 

- R £sin 9 cos 6 c cos|^ c - </: 

~ coshcos(? c J , 

(12) 



- c 

• 0 = 

R sin 6 sin - ./-) 





Recalling that 

B r f + Bg 0 + } = B -- W 

where the potential V is given by the well known function 


v = 


L 


(t) 1 ('O (g n m cos# + h™ sin •;) , 


(14) 


we obtain the components of B : 


dV 
d r 


j_iy 

r d t ‘ 

1 dV_ 

r sin ' rV 


(15) 


for increasing r, andc/ . The coordinate system used is illustrated in Figure 5b. The field com- 
ponents along the unit vectors u v ; - ? , u N T - o and u E / are 


B f : vertically downward 


Bo : horizontally north 


B E = B^ : horizontally east 


(16) 
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The field calculations are performed with Mcllwain’s computer program MAGNET, using the 
99-term Hendricks and Cain model for 1960.0 (GSFC-9/1965) updated to 1965.0. From the value 
of the field strength and the geometric considerations above, an accurate evaluation of Equation 8 
is possible. 

RESULTS 


The results of our computations are given in Figures 6, 7, 8, 9, 10, and 11 in terms of par- 
ametric angles vs. time. Figures 12 and 13 give the spot-to-station distance with time. For a fixed 
altitude H, the look-elevation e is strictly a function of the line-of-sight distance D: 


e - arccos ( 1 - H^D 2 ) 1/ ' 2 


where D = |c| . 

Figure 14 is a plot of e vs. D for H = 100 km. It may be used for quick determination of the 
approximate elevation, if the distance is known. Predictably, whenever e is greater than 88°, the 
look-azimuth 0 becomes meaningless as a parameter for tracking. This condition exists when the 
intersect trace passes very close to, or over, a station. The effect on 0 of a near zenith passage 
of the spots is an abrupt jump in quadrants; it poses problems in adjusting instruments that cannot 
be moved through their local zenith. Stations No. 6 and 7 are affected most, whereas Stations No. 

1, 2, 8, and 9, located well to the north or to the south along the extended projections of the main 
intersect segment, experience minimum variations in 0. The consequence of these conditions is 
reflected, for the combined data of e and 0 from all Stations, in Figures 15 and 16, which depict 
respectively "optimum" and "critical" look-angle areas on the ground. 

In the context of this study, artificial auroras are spots generated by a rocket-borne accelerator 
that emits electron beam pulses of constant energy at intervals of three seconds. The electrons, 
interacting with the constituents of the upper atmosphere, are expected to produce the spots at an 
altitude which depends largely on their kinetic energy (Berger et al., 1968); the more energetic 
particles penetrating deeper into the atmosphere, the less energetic interacting at a greater height, 
(Figure 17). The sequential appearance of the spots in time should follow closely the intersect 
trace of that altitude level. 

Corresponding to every rocket-generated intersect trace there exists an area on the ground 
where the auroras will be brightest for most of the duration of the flight above 150 kilometers, the 
luminance at a given location being proportional to a factor k = l/sin a. Figure 18 shows the 
value of the luminosity factor over the range 0.5° < a <90°. 

Using the information available from Figures 10 and 11, the optimum o'-area was determined 
on Figure 19 for k > 2 and k > 3. The apparent brightness of the spots may be of interest to some 
experimenters when selecting the site for a camera station. 
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From the evaluation of Figures 15, 16, and 19 the area of "Optimum Parametric Angle Condi- 
tions" emerges and the "Critical Range" develops, as shown in Figure 20. In the "Optimum Area”, 
I corresponds to a luminosity factor of k > 2 and n to k >3. 

Finally, a station-centered polar look-angle plot is added as in Figure 21. It represents the 
intersect trace for each station plotted in e vs. \p. 
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Appendix 


Field Line Tracing 


The coordinate system used in the field line tracing process is illustrated in Figure 5b. The 
field components along the unit vectors u v , u N , and u E are: 

B v : vertical down 


B n : horizontal north 


B E : horizontal east 


The equations defining an infinitestimal portion of a line of force are: 


ss 

B 




cos K&4> 


and SS is the element of arc, with components - Sr, rSA. f and r cos \S<j> along the unit vectors. The 
coordinates of a generic point of a field line originating at r s , \ s , and <t> s are then given by: 


n B 


r - r + <j 


k - A . - a 


IX ss * ■ 

i = 1 

y -1 1 Bn i 


</>„ - <t>. - CT 


Op 

r 1 l E i „ 

/ r. cos k . B. 5S i 

/ -i 1 1 1 


where the sign factor o- is ±1 according to whether one proceeds tracing in the same direction as 
the field vector or the opposite. If the origin is in the magnetic southern hemisphere (B v < 0), it 
is necessary to set a < 0 in order to follow the descending line to the nearest intersect. 
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A four point "Adams integration formula" for numerical solutions of differential equations 
is used in the tracing procedure. Its repeated application at equidistant intervals approximates 
the field line. The accuracy of this fit depends on the integration step size. 
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FIGURES AND CAPTIONS 


Figure 

1 Range of magnetic shells crossed by Aerobee trajectory. 

2 Map of stations with trajectory projection and 100-km intersect trace. 

3 Geometry and relationships employed in evaluating the azimuth ^ina spherical 
coordinate system. 

4 Geometry of vectors and angles used in evaluating the elevation e in a spherical 
coordinate system. 

5 Geometry of vectors and angles used in calculating the line angle a. 

6 Variation of look-elevation with time: Stations 1 to 9. 

7 Variation of look-elevation with time: Stations 10 to 19. 

8 Variation of look-azimuth with time : Stations 1 to 9. 

9 Variation of look-azimuth with time: Stations 10 to 19. 

10 Variation of line-angle with time: Stations 1 to 9. 

11 Variation of line-angle with time: Stations 10 to 19. 

12 Variation of distance with time: Stations 1 to 9. 

13 Variation of distance with time: Stations 10 to 19. 

14 Distance dependence of look-elevation for points on the 100-km. altitude shell. 

15 Area of optimum look-angle conditions. 

16 Critical e and \p area. 

17 Approximate interaction altitude of electron beam. 

18 Values of luminosity factor over the range: .5° < a < 90°. 

19 Optimum a area. 

20 Area of optimum parametric angle conditions and critical range. 

21 Station- centered polar look-angle plot. 
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Figure 2— Map of stations with trajectory projection and 100-km intersect trace 









Figure 6— Variation of look-elevation with time: 


Stations 1 to 9. 
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Figure 7— Variation of look -elevation with time: Stations 10 to 19. 
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Figure 9-VariaHon of look -azimuth with time: Stations 10 to 19. 
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Figure 11 —Variation of line-angle with time: Stations 10 to 19. 
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Figure 16— Critical e and area. 
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